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Extreme value statistics of 2d Gaussian Free Field: effect of finite domains 


Xiangyu Cao, Alberto Rosso, and Raoul Santachiara 
LPTMS, CNRS (UMR 8626), Universite Paris-Saclay, 91405 Orsay, Franc^ 

We study minima statistics of the 2d Gaussian Free Field on circles in the unit disk with Dirichlet 
boundary condition. Free energy distributions of the associated Random Energy models are exactly 
calculated in the high temperature phase, and shown to satisfy the duality property, which enables us 
to predict the minima distribution by assuming the freezing scenario. Numerical tests are provided. 
Related questions concerning the GFF on a sphere are also considered. 


I. INTRODUCTION 


The minima statistics of the 2D Gaussian Free Field (GFF) is a fundamental problem relevant in different contexts, 
ranging from multi-fractal behaviours in quantum disordered systems EHi , rare events in climate science Q to zeros 
of Riemann zeta function Q . While rigorous results exist for the leading behaviour of extrema and near extrema of 
GFF, c.f. a and references therein, less is known about the full extremum distribution. This was calculated only in 
two cases: the GFF on a circle by Fyodorov and Bouchaud Q (FB), and on an interval 3- Fi both cases, the authors 
studied the associated Random Energy Model ^ O '''^Fh logarithmically correlated energy landscape. They showed 
that the minima distribution can be obtained by a freezing seenario describing the low temperature phase of those 
Random Energy Models [HI- Moreover the authors of S pointed out the key role of a duality property, reminiscent 
of well-known dualities in Gonformal Field Theory (l^. [l3|. 

In general, the GFF 4>^{x) on a 2D domain E is determined by its Green function 4 >^{x)4)y.(]j) = G^{x,y) ~ 
—2\x\\x — y\. It has both a short distance (ultraviolet) and a long distance (infra red) divergence. While the first 
is at the origin of the freezing scenario, the second must be cured in some way. For instance, one mw take E to 
be the unit disk and specify the Dirichlet boundary condition. From this point of view, the results in [7| should be 
interpreted as the limit of vanishing circle radius. In this work, we study GFF on circles of any radius, where the 
minimum distribution is expected to be affected by the presence of boundary, c.f. sect. 6.2.1. Alternatively, E 
can be a compact surface, such as a sphere. One of our main motivations is to see whether freezing and duality are 
robust in these different situations. 

When E is a compact surface, the Green function ^ is the solution kernel to the Poisson equation. For the sphere 
E = §, with the standard metric, we have (see e.g., 


Gs{z, w) = — In 


\z — wl 


(1 + |^|2)(I + |^|2) 


- 1 


( 1 ) 


where we identify S with C U {oo} by stereographic projection. On a latitude {j^j 


= r} and noting (j)o{9) := (/)s(re*®), 


4>o{9)(t)o{9') = -21n|e‘® - e‘® | -I- a. 


( 2 ) 


Here a is a zero mode, which contributes a convolution with an independent Gaussian in terms of minima distribution. 
We shall set a to zero. Then the minima distribution of the GFF on a latitude on S does not depend on r, but always 
coincides with the FB case. 

As we will see, this is no longer true for the unit disk ID), the Green function on ID with Dirichlet boundary condition 


is [I1 ]Gd( 2 :,ui) = -2 In 


Z — W 

1 — zw 


. On a centred circle of radius q € (0,1), noting (f(0) := ^^(y^e'®) 


(t){e)(j){e') = -2 In 


I - z 
1 — qz 


z:= 


(3) 


where a' is the zero mode that we will set to zero as before. In terms of Fourier modes, 

(j,(0) = y/)I^exp(i6>fc)Nfc, pik := 7^(1 - g'^'), k 0 


( 4 ) 
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Here are complex numbers whose real and imaginary parts are i.i.d standard normal random variables, and JRz is 
the real part of 2 . When q ^ 0 the model eq. ([3]) reduces to the FB case Note that the model defined by eq. 
(|4|) makes perfect sense also for g G [—1,0). For this range of g, eq. ([3]) can be obtained from the Green function 
Gp{z,w) = Gs{z,w) — Gs{z,w) (where w = —w~^ is the antipodal of w), by taking | 2 ;| = jwj = ^y—q. We can then 
interpret the results for g < 0 as the minimum of := {4>s{z) — 4)%{z))/'/2 on a latitude of S. In particular, for 

g = — 1, we are considering the extrema of [4>o{0) — if)Q{9 + tt))/\/ 2, (j)o being that of the FB model. 

In the sect. m we will solve the model defined by eq. ©; sect. imi is devoted to numerical tests of the solution. 


II. DISTRIBUTION OF MINIMA 


Let us begin by clarify the term “minima of GFF”. Indeed, one dehnes first the minima of finite systems ijm '■= 
Here fi = and (j) is defined by eq. o where the sum is cut to jfcj < We are interested 

in the rescaled minima y = limM->.oo > where the scaling behaviours cm = —21nM + |lnlnM + 0(1) and 

&M = 1 + O were predicted in [n|. 

To calculate the distribution of y, we study the Random Energy Model defined by the partition function Zm = 
/3 = T~^ is the inverse temperature. Its positive moments (or replica averages) (n = 2, 3 ...) can be 
calculated by applying Wick theorem and using eq. For each n, for /3 small enough, the sum can be approximated 
by a Goulomb gas integral as M — 00 : 


where 


'P'n _L. 
— 



n 




Zri. ' — 


ntn 

i<j 


Zi - Zj 


1 - qZiZj 


-27 


7:= P - 


(5) 


( 6 ) 


This integral converges for /3 < /3c,n = n and this is precisely the condition for eq. m to hold, c.f. 0 sect 
3.1. The temperature /3c,rt’s, called pre-freezing temperatures, are precursors of a true transition separating a high-T 
phase and a glassy phase. To see this heuristically, one can compare the number of energy levels, M, and the mean 
Boltzmann weight of each level, oc [s^ A transition occurs at Pc = Pc,i = 1: 

- When /3 < 1, M > , Zm is entropy dominated. So we expect that free energy distribution in the 

thermodynamic limit can be obtained by analytically continuing the Goulomb-gas integral ([SI), which we do in 
sect. MM 

- When /3>1, > M, Zm is dominated by a few (^ 0(1)) low energy levels and the leading term (oc InM) 

of the entropy vanishes. This non-analyticity extends to the full free energy distribution, and invalidates any 
analytic expansion method. The freezing scenario allows to extend non-analytically ,5 < 1 results to ,5 > 1, and 
to predict the distribution of y at /3 —>■ 00 , see sect. ME 


A. Analytic continuation of moments (when /3 < 1) 

The analytical expansion Zn will be done using Jack polynomials [l^ (following conventions of [13). Denoted by 
Jx{z; a), they depend on variables z = (zi,..., Zn), a parameter a and a partition A. The latter is essentially a Young 
diagram, i.e., a set of unit square boxes with integer coordinates, c.f. fig 0 

A = {(a;,?/) : 0 < a: < Ay,0 < y </(A)} (7) 

where Aq > Ai > ... > Ai(A)_i > 0 {1{X) is its length). We denote by |A| := f^e total number of boxes in 

A, i.e., its size. 

An important property of the Jack polynomials is the Cauchy identity [isj] 

n 

(1 - qZiWj)-^ =^q^^''Jx{z-,a)Jx{w,a)jf^{a), 

i,j=l X 


( 8 ) 
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where the RHS sums over all the partitions. The value of j\{a) will turn out irrelevant. Setting a = — 7 , Zi = 
and Wi = , we can use eq. dH) to rewrite the denominator of integrand of eq. dH]) as 


n 


1 

1 - qzizj 


-27 


Jx{z;a)Jx{z;a) 

, 3x{a) 


where z = {zi,..., z„). We shall then need the inner product norm 


(9) 


where A := ]\i^j{zi 


TT 2 a J\iz; a)J\{z\ a) 

Jo l\ 2n ' ' c„(a)jA(a) 


n 9 n{x,y;a), 
{x,y)e\ 


Cn(a) 


TT ddi I A |2a _ r(l + na) 

Jo " r(l + a)- 


gn{x,y;a) 


X + a(ri — y) 
a; + 1 + a(n — y — \) 


Zj), and Cn is the Dyson integral [ 2 ^. Combining dH) and (fTUl) . we obtain 


( 10 ) 

( 11 ) 

( 12 ) 


^n(/3,<?) = (1 - < 7 ) "^c„(-7)^g*^l 5„(x,?/;-7) 

A (x,y)^\ 


(13) 


This equation makes sense for n and /3 generic, which allows us to extend the definition Z„(/3, q) beyond the domain 
where corresponding Coulomb-gas integrals exist. 

Now we define the fluctuation part / of the free energy Fm = —InZjvf as 


/:= 


-/3 Mnz, z := Ze := Me 


X^(l-9)' 


r(i-7) 


(14) 


So / and Fm differ by a shift, containing the M-dependent part, so that the distribution of / has a well-defined limit 
as M —>• 00 . Indeed, writing t = —n/3, one verifies that the equations dS]), and (ITT)) imply (at M = 00 ) 


s(t,/3,g) := 

A 


n 

(x.i/)GA 


expt/ = T{l + tf3)s{t,l3,q), 
xP + yP~^ + t 
[x + l)f3 + {y + l)/3“^ -I -1 ' 


(15) 

(16) 


Note that when g = 0, s = 1 and expt/ = r(l -|- t/?), coinciding with the FB solution. The series s(t,/3,q) encodes 
the effects of the finite domain, and is the main result of this work (supplemented by an efhciently calculable matrix 
product rewriting, c.f. Our approach relying on formal analytical continuation of eq. m and m is non-rigorous 
and the result eq. (fT3)) is a conjecture. Its validity is well supported by careful numerical tests (sect. IIIip . Moreover, 
in m we check that it reproduces correctly a non-trivial coefficient in the hi gh-T expansion. To prove rigorously eq. 
(ESI we should compute the high-T expansion to all order, similarly as in [^1^. but this project is left to the future. 


B. Freezing and duality (when /3 > 1) 


Eq. (1151) predicts free energy fluctuations for /3 < 1. For /3 > 1, we shall assume the freezing scenario. It appeared 
first as front velocity selection criteria in Fisher-Kolmogorov-Petrovsky-Piscounov type equations [1^, HJ. Later it 
emerged in the physics of disordered systems: first these defined on trees 1^. 1^ . and then Euclidean-space ones 
[l|, 0, [li|. Recently, it is proved by Madaule et. al. for a large class of models [^. To state it in our context, consider 
the random variable yp defined by yp '■= f — (3~^g, where g is a Gumbel variable independent of /. Equivalently, 

exp(%) = exp(t/)F (^1 + ^ (17) 

By eq. [151 the above has a well-defined limit as /3 1. Now the freezing scenario claims it becomes /3-independent 

for /3 > 1: 


exp(%) = eycpityp) 


p/'i 


(18) 
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Figure 1: Analytical prediction for minima distribution. Left: The mean is shifted to zero. Right: The distributions are 
rescaled to have zero mean and unit variance. 


In particular, taking the /3 —>■ oo limit, we have yp=i = yp^oo = /(/? —t oo) = y, giving the minimum distribution. Its 
p.d.f., calculated by inverse-Fourier-transforming eq. m is plotted in fig. ©• Remark that s(t, /3 = 1, g) have poles 
only at t = —2, —3,..., (eq. (ITB)) 1. so the rightmost pole of exp(ty) is at t = —1 and of order 2, (eq. (fTTll l. By inverse 
Fourier transform, this implies Pdf(y —>■ —oo) ^ for any q, confirming the universality of the Carpentier-Le 

Doussal tail [HI- 

We turn now to the duality. Authors of Q observed for the FB model and the interval model that, if one continues 
analytically expityp) to /3 > 1, the unique result is self-dual, i.e., invariant under j3 :—>■ /3“^. We stress that the 
self-dual solution is physically wrong for the phase ft > 1, where the non-analytical continuation eq. m holds. An 
important result of this work is that, exp^typ) enjoys the same duality property, unaffected by finite domain effects. 
It follows from eq. (fT5|) and the self-duality of s(t,/3,g) of eq. (fTHll : 

s{t,(3,q) = s {t,P~'^,q) . 

In fact, the term for A and for its transpose A' = {{y,x) : (^y) G A} are related by /3 O fi~^. The co-presence 
of freezing and duality property supports the conjecture of [Sj , which remains nonetheless intriguing from a general 
theoretical viewpoint. 


III. NUMERICAL STUDY 

Numerical simulations of our models follow the same protocol as in Q, sect. 6.1. The finite models defined at 
the beginning of sect. |lT]are simulated using Fast Fourier Transform. We will focus on: (i) the validity of analytical 
continuation leading to eq. (US and (ii) the freezing scenario by measuring the free energy variance as a function 
of temperature. 

We first study the free energy distribution in the ft < 1 phase, predicted by eq. (HID. To compare analytic 
and numerical values, the mean value should be shifted away from both. For numerical data, we measure the 
empirical mean value and subtract it from each sample: /' := fi — (X]J=i ^.nd calculate the Fourier transform 

^"^1 exp(i//t)/n of the sample distribution. On the analytical side, we calculate exp(i/t) exp(—i/t). The sum 
s{t,ft,q) in eq. dTSl) is calculated with the matrix product eq. (IA 2 I) . and the first moment / using eqs. (IA5I) and 
(IA 6 D . In fig. ([2]), we plot real and imaginary parts of the Fourier transform, so t becomes an invisible parameter. 
This amounts to rescaling the distributions to have variance 1, so that we concentrate here on the detailed shape of 
the distribution, e.y., its asymmetry. The results validate the prediction of the free energy distribution in the high-T 
phase. 

Now we test the freezing scenario. We follow a strategy of Q (sect. 6.4) and look at the variance of free energy 
as a function of ft. The ft < 1 part supplements the precedent numerics, while the ft > 1 part tests the freezing 
scenario proper. The analytical prediction is calculated from the matrix product eqs. (IA5I) . (IA 6 I) and (IA7I) for ft <1, 
and continued according to eq. ITSl) to ,5 > 1 . The results are shown in fig. Although convergence in the ,5 > 1 
phase is notoriously slow, we observe tantamount evidence of freezing for all values of q. The validity of the predicted 
free energy distribution in /? < 1 phase, plus the freezing scenario, entail the correctness of the prediction minima 
statistics. 
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Figure 2: Analytic vs. Numerical calculation of Fourier transform of free energy exp(it/), t > 0. Besides what is plotted, 
several other combinations (q, /3) are also checked. Left: In the high temperature phase, here ft = 0.5, the analytic-numerical 
agreement is excellent for all values of q from —1 up to 0.9. Right: At the transition, finite size corrections amplify, especially 
for q 1. 



Figure 3: The variance of with the FB {q = 0) analytic values subtracted. The analytic-numerics agreement is excellent in 
the high-T phase and degrades when approaching the transition. In the low temperature phase, the plateau indicating freezing 
is observed for all numerical data. 


IV. CONCLUSION 

In this work, we mainly focused on the GFF on the unit disk with Dirichlet boundary condition. We proposed 
the minima distribution of the GFF on a centred circle of radius ^/q for general q € (0,1) (fig. ([1])). This was done 
by first proposing the exact free energy distribution of the Random Energy Model for /3 < 1, tested with excellent 
agreement against numerics, fig. (ED. and then applying the freezing scenario, eq. CHD. Furthermore, as explained in 
sect. IIIBI the self-duality enjoyed by our solution. Interpretation of results when q G [—1,0) and related problems 
on the sphere is also discussed (sect. |T|. The results presented here are the first that probe the effects of boundary 
conditions or finite volume. We show that these latter do affect the minimum distribution, but keep intact duality 
and freezing features. Their co-presence is observed for exp(tj/^) of all solved 2d GFF models, as well as for velocity 
moments in decaying Burgers turbulence [2^. It is interesting to investigate its gen erality, origin and relation with 
the duality appearing in Random Matrix Theory [s^ or Gonformal Field Theory [1^ . 

Acknowledgements. We are grateful towards P. Le Doussal for a careful reading of the manuscript. We thank E. 
Bogomolny, A. De Luca, O. Giraud and G. Texier for useful discussions. 


Appendix A: Matrix product form 

The inhnite sum s(t,ft,q) of eq. (1161) is absolutely convergent for jqj < 1 and away from poles, but inefficient for 
practical use. Moreover, for q = —1, s is not absolutely convergent, and the naive partial sums X]|a|< 7 V oscillate 
> 10% even when N ^ 80. Fortunately, both problems can be resolved by a matrix product form of s(t,/3,q), which 
we derive here. 

First, we consider partitions as lattice walks ft,(r),T G Z, where hlr) is the number of squares on the diagonal 
X — y = T (see fig. (|3D for explanation). Observe that the product over boxes in eq. (fTO telescopes along diagonals 
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h 


In general, /i(r) satisfies: 

a. /i(r) — h{T + 1) e {0,1}, r > 0; 

b. /i(r) — h{T — 1) G {0,1}, T < 0; 

c. limT-»±oo /i(t) = 0. 


Figure 4: Left-. An example illustrating the mapping between partitions A and paths Right-. Conditions h should satisfy 

to correspond to a partition. 


and becomes a product over t. We have indeed 


s(t,/3,g) 


E n 


g^(^)(|T|/3"sn(^) +t) 
|r|/3"gn(^) + (/3 + /3-1) /i(r) + t 


(Al) 


The sum is over all walks /i(t) satisfying conditions a, b and c of fig. (|H). The RHS resembles a path integral (or the 
partition function of a directed polymer), and can be written as a product of (transfer) matrices. To avoid infinite 
vector spaces, consider the truncation of ED to paths h{T),T = —I,... ,l such that h{±l) = 0 (j.e., we sum over 
all Young diagrams contained in the square {x < l,y < ?}). Using the auxiliary space spanned by |h), h = 0,1, 2 ... /, 
we have 


3 W = {0\UDi_iU ... DiU\C\LEi... LEi_xL\Q) 


1-1 


U = I + J2\h){h + llD,{/3) = Y. 


q^\h){h\ 


h^O 


l + h(/3 + ^-i)(j^ + t) 


-1 ■ 


L = u\Ej= d} (/ 3 - 1 ) ,C = Do = Eo. 


(A2) 

(A3) 

(A4) 


Remark that Dj,Ej and C are diagonal (they generate factors in eq. (lAll) . j = |t|), while U and L are nearly so 
(they implement conditions a and b in fig. (|T]). This enables the sum s{l) over 0(4*) partitions to be calculated in 
0{P) time and 0(/) space, so as to achieve in practice the convergence s, I /' 00 . Using I ^ 10^, we observed 

4 decimal precisions, even for q = —1. Below, s will be understood as for I sufficiently big. 

To calculate moments of /, we need t-derivatives of eq. (IA2I1 at t = 0. In fact, from eq. CSl) we have 


Sfc:= ^ K = Sfc + ^"(lnr(a:))i"2i, P<1. (A5) 

For example, /^ = si — /37 b, = S 2 + /3^7r^/6, and so on. Differentiating s is done by Leibniz rule, and the result 

is written again as matrix products. For the first moment, one has 

Si = {0\UdI°_\u ...D^°'>U\C’-^^LE[°'> . ..LeI°_\l\0) (A6) 


where = Dj\t=o, similarly for Ej^\ and = {dtC)\^^Q. 

For the fe-th moment, one need an auxiliary space of (/ + l)fc-dimension (to bookkeep the positions hit by derivatives), 
spanned by |/i, d) := \h) ® |c?), h = 0,... ,l,d = 0,... ,k — 1. The general expression is cumbersome, and we will give 
that of = 2: 

S 2 +sl = {0,0 \UdI]}^U ... d\^^U\C‘-‘^^LE[^'’ ... LeI]_\l\0,0) (A7) 

dW + |0)(l|,U = Dt(/3^/3-i) (A8) 

=2(atC)Uo®(|0)(l| + |l)(0|)+ (92c)2^^®|0)(0| (A9) 
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Appendix B: High temperature expansion 


As an another check of eq. (HSJ, we compute the high temperature expansion of the free energy variance at 
order/3^. As is shown in @ , (eq. (C.5)),We have= (1 — (eq. S]), which implies 

F = (Li2(q2) _ 2Li2(g) + y) /3" + Li™(x) := E (Bl) 

^ fc >0 ^ 

Now we derive this from eq. m- The r(l + tf3) part gives the contribution The sum eq. (flBl) is developed at 

t = /3 = 0 as: 9 '^' n(a;,y)eA ^ + t{mif3 + ...)+ where ... are higher order 

terms. 

Let us begin by mi. For any non-empty partition, the box (a; = 0,?/ = 0) has already nominator t/3; so the other 
boxes should all be evaluated at t = /3 = 0 when contributing to mi. So, only A’s with one column contribute; 
otherwise the box (x = 1 , 2 / = 0) would have a nominator (/3^ + t/3) that vanishes. Therefore 


fe-i 


,,-n ;c>0 


(B2) 


fe>0 y—0 


Similarly only one-column and two-column partitions contribute to m 2 . A two-column partition contributes by the 
product of all boxes with y > 0 evaluated at /? = t = 0, yielding 


k-l 


i-i 


m 2 


, 2 - E '?'^'n^n^-^Lii(gf+iLi2(F 


k>l>0 


(B3) 


y=o 


Finally, contributions from each one column partition is the sum over where the second t/3 term comes from: 


m 2 


■i-Z-T 


fe >0 


k-l 


-' + 1 : ^ 


^ vy 2/ +1 


= -Li 2 (g) 


(B4) 


Combining eqs. IB2l to [B4l one has 2m2 — mf = 2(m2,i + m 2 , 2 ) — mf = Li 2 (g^) — 2Li2(g) as desired. 
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